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Abstract 

We study the effects of the finite number of experimental data on the computation of a 
generalized fluctuation-dissipation relation around a nonequilibrium steady state of a Brow- 
nian particle in a toroidal optical trap. We show that the finite sampling has two different 
effects, which can give rise to a poor estimate of the linear response function. The first 
concerns the accessibility of the generalized fluctuation-dissipation relation due to the finite 
number of actual perturbations imposed to the control parameter. The second concerns the 
propagation of the error made at the initial sampling of the external perturbation of the 
system. This can be highly enhanced by introducing an estimator which corrects the error 
of the initial sampled condition. When these two effects are taken into account in the data 
analysis, the generalized fluctuation-dissipation relation is verified experimentally. 

1 Introduction 

Current theoretical developments in nonequilibrium statistical mechanics have led to significant 
progress in the study of systems around states far from thermal equilibrium. Systems in nonequi- 
librium steady states (NESS) are the simplest examples because the dynamics of their degrees of 
freedom x under fixed control parameters A can be statistically described by time- independent 
probability densities po{x,X). NESS naturally occur in mesoscopic systems such as colloidal 
particles dragged by optical tweezeres, Brownian ratches and molecular motors because of the 
presence of nonconservative or time-dependent forces [1]. At these lengthscales fluctuations are 
important so it is essential to establish a quantitative link between the statistical properties 
of the NESS fluctuations and the response of the system to external perturbations. Around 
thermal equilibrium this link is provided by the fluctuation-dissipation theorem [5]. 

The generahzation of the fluctuation-dissipation theorem around NESS for systems with 
Markovian dynamics has been achieved in recent years from different theoretical approaches 
[3l[l[5l|6l[3IHl[9l[10llIU[12]. The different generalized formulations of the fluctuation-dissipation 
theorem link correlation functions of the fluctuations of the observable of interest 0{x) in the 
unperturbed NESS with the linear response function of 0{x) due to a small external time- 
dependent perturbation around the NESS. The observables involved in such relations are not 
unique but they are equivalent in the sense that they lead to the same values of the linear 
response function. These theoretical relations may be be useful in experiments and simulations 
to know the linear response of the system around NESS. Indeed the response can be obtaines 
from measurements entirely done at the unperturbed NESS of the system of interest without 
any need to perform the actual perturbation. Nevertheless, the theoretical equivalence of the 
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different observables involved in those relations does not translate into equivalent experimental 
accessibility: e.g. strongly fluctuating observables such as instantaneous velocities may lead to 
large statistical errors in the measurements [13]. Besides, NESS quantities themselves such as 
local mean velocities, joint stationary densities and the stochastic entropy are not in general 
as easily measurable as dynamical observables directly related to the degrees of freedom [H] . 
Hence, before implementing the different fluctuation-response formulae in real situations it is 
important to test its experimental validity under very well controlled conditions and to assess the 
influence of finite data analysis. The experimental test of some fluctuation-dissipation relations 
has been recently done in Refs. [131 [TH [121 [E] for colloidal particles in toroidal optical traps. 

In the present paper we discuss the effects of the finite number of experimental data on 
the determination of the linear response function around a NESS for a micron-sized system 
with Markovian dynamics: a Brownian particle in a toroidal optical trap. For this purpose we 
perform the respective data analysis on the measurements reported in [Ml [T6|. In Sect. [2] we 
briefly describe a generalized fluctuation-dissipation relation that has been derived for Markovian 
dynamics around a NESS exploiting the properties of the stationary density po{x, A). In Sect. [3] 
we recall the main features of a previous experiment that we use in the present paper for the data 
analysis. In Sect.jJlwe discuss the two different kinds of finite-sampling effects that can appear in 
the computation of the differents terms involved in the fluctuation-dissipation relation. We show 
that the generalized fluctuation-dissipation relation is verified experimentally when performing 
a careful data analysis, which takes into account these effects. Finally we present the conclusion. 



2 Hatano-Sasa relation and fluctuation-dissipation around 
NESS 

The Hatano-Sasa relation provides a general identity for the transitions between either equilib- 
rium or nonequilibrium steady states of Markovian systems [T7] . In the following we will focus 
on a Langevin system with Markovian dynamics described by a steady state probability density 
Po{x, A). When the system is subjected to a time-dependent variation of the control parameter 
\{t) between an initial time U and a final time i/, the Hatano-Sasa identity reads 

where (j){x,X) = —\i\pq{x,X) and the average (...) is performed over an infinite number of 
realizations of a prescribed time-dependent protocol \{t). From Eq. ([T]), Prost et al. has directly 
derived a generalized fluctuation-dissipation relation that holds in the linear response regime 
around a NESS [7] 

d / d(j){xt,Xss) d(l}{xs,\ss)\ 

''-^'"^^^d^V^A; dx;—/,- 

In Eq. ([2]) Ra-y{t — s) = S{0{xt))h/Shs\h=o is the linear response function of the observable 
0{x) = d(j){x,\ss)ld\a due to a small external time-dependent perturbation hg ~ A(s) — Ass 
around Ass = A(ti) fixing an initial NESS at time ti. The averages {...)h and (. . .)o are 
performed over the perturbed and unperturbed processes, respectively. 

In experiments the formal average involved in Eq. ([T|) is not perfectly computed because of 
the finite number of independent realizations of \{t). Hence Eq. ([1]) allows one to estimate the 
experimental precision of ^ computed from a given number of experimental data provided that 
one can measure the observable d(j){x, Xss)/dXa- In the next section we tackle this problem for 
the experimental trajectories of a colloidal particle in a toroidal optical trap. 
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Figure 1: (a) Sketch of a Brownian particle in a toroidal optical trap subjected to a nonconserva- 
tivc force / and undergoing a periodic potential U{9).{h) Experimental profile of the observable 
0{9) defined in Eq. (g]) computed using the NESS densities poi9,A + SA) and poi9,A - 6 A) 
around po{9,A) at fixed F. 

3 Brownian particle in a toroidal optical trap 

The Brownian motion of a colloidal particle in a toroidal optical trap has become an experimental 
model to study the generalization of the fluctuation-dissipation theorem around a NESS [13 1 [M l 
[TSl [T6j . This is because it is a system with a single translational degree of freedom where one can 
easily tune its relevant control parameters. Our experiment has been already described in detail 
in |14[ I16j so here we only explain it briefly. The Brownian motion of a spherical silica particle 
(radius r = 1 pm.) immersed in water is confined on a thin torus of major radius a = 4.12 /zm by 
a tightly focused laser beam rotating at 200 Hz. The water reservoir acts as a thermal bath at 
fixed temperature (T = 20 ± 0.5°C) providing thermal fluctuations to the particle. The viscous 
drag coefRcient at this temperature is 7 = 1.89 x 10^^ kg s~^. The rotation frequency of the laser 
is so high that it is not able to trap continuously the particle in the focus because the viscous 
drag force of the surrounding water quickly exceeds the optical trapping force. Consequently, 
at each rotation the beam only kicks the particle a small distance along the circle of radius a 
exerting a nonconservative force / = 66 fN on it in the direction of the rotation. Thus, the 
particle motion is effectively confined on a circle: the angular position 9 of its barycenter is the 
only relevant degree of freedom. In addition, a static light intensity profile (amplitude about 
5% of the total laser intensity [Ml [1^ ) is created along the circle acting as a periodic potential 
U{9) = U{9 + 2tt) of amplitude 68.8kBT. Figure [Ija) depicts this experimental configuration. 
We track the 2D particle position by video microscopy in order to measure the time evolution of 
9. Thus for the experimentally accessible length and time scales the dynamics of 9t is modeled 
by the first-order Langevin equation [T31 HH [13 [HI fT5] 



where H{9) ^ U{9)/{ja) with amplitude A = raa^{H{9)} ^ mMsT / {"fo^), F = //(7a) and 
^ is a white noise process of zero mean and covariance {£,tS,s) = '^[ksT / {"fa?)]&{t ~ s). Under 
these fixed conditions the dynamics of 9t settles into a NESS whose probability density function 
pa{9,A) is plotted in Fig[Ub) (solid black line). 
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4 Generalized fluctuation-dissipation relation 



In the following analysis we take a; = 6' as the single degree of freedom and A = A as the 
main control parameter of the system. In this case the observable of interest involved in the 
fluctuation-dissipation relation ^ is 

OW ^ (4, 

The experimental profile of 0(9), computed as a discrete three-point derivative of — lnpo(^5^) 
at two different NESS around A, is shown in Fig[ljb). 

We focus on the response of the system after applying a small Heaviside perturbation SA 
to A at time ti: A ^ A + 6 A. In the experiment this dynamical procedure is done by sud- 
denly switching the laser power modulation as explained in |14| I16j . This procedure yields the 
integrated response function, defined as 

x{t-t.) = j\{t-s)ds = ^^^MlA, (5) 

where the average (. . .)sa rnust be performed over the perturbed process at time t. Then the 
integrated version of the generalized fluctuation-dissipation relation ([2]) is in this case 

X{t - U) = {O{et)0{9u))o - {0{9t)0{dt))o. (6) 

We now study the effects of a finite number of realizations of the perturbation SA and the finite 
number of trajectories used to compute the averages (. . .)sa and (. . .)o in Eqs. ([5]) and © . 

4.1 Statistical error 

As discussed in Sect. [2l the Hatano-Sasa relation ([T]) can be used to estimate the error of the 
experimental computation of Eq. (jG]) when performing N < oo independent realizations of SA 
around the NESS. In the case of the dynamical process defined by the Heaviside perturbation 
A ^ A + SA at time U, as done in the experiment, Eq. ([T]) reads 

{exp[-SAO{eu)])^l. (7) 

Eq. ^ only depends on the initial values 0{6ti) when the system is still in NESS. Therefore, 
for a finite number of trajectories N we introduce an estimator of the error of Eq. ([7]) 
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(8) 



where Oj{9t-) is the j-th sampled NESS initial condition. Fig. [2^a) shows the behavior of the 
error A (TV) computed for SA = 0.05 A (the value realized in the dynamical experiment) using 
N experimental values of 0{9) drawn from the NESS distribution. For small N < 100 the error 
is non-ncgliglible, A(iV) > 4%, and this must be taken into account in the final accuracy of 
the generalized fluctuation-dissipation relation when comparing the left with the right-hand side 
of Eq. (jni) using the experimcnal data. Then as N increases A(A^) quickly converges to 0: for 
N > 500 the precision of Eq. ([7]) found in the experiment is better than 1%. 
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Figure 2: (a) Estimate of the error of the Hatano-Sasa relation ([7]) for a finite number N of 
realizations of 5 A. (b) Estimate of the error of {O{6tj)o using Eq. ^TU^ for M NESS data. 



4.2 Effect of the initial sampled condition 

According to the Heaviside procedure for SA, in Eqs. (O and ^ the initial condition 9{ti) for 
the perturbed process is sampled from the NESS density pq{6,A). Then the integreated linear 
response function must formally satisfy the initial condition 

X(0) = = 0, (9) 

where the last equality is due to the normalization of po{9ti,A). We are interested in the effect 
of a finite number M < oo of initial values 9{ti) drawn from the inital NESS on the estimate of 
xit — ti). It should be noted that in practice a small M may significantly affect the computation 
of {O{9ti))o in Eq. ^ because most of the positive values of 0(9) he in the region where po{9, A) 
is rarely sampled, as shown in Fig. [IJb). In Fig. ^b) we plot some values of the finite average 
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where Oj{9t-) is the j-th initial condition at NESS. As expected, for small M, {0{9ti))M < 
due to the fact that one samples mostly the negative values around the maximum of po{9,A). 
The convergence to the theoretical value (O(6'tj)o = is very slow: as M increases {0{9t-))M 
becomes very sensitive to M and large positive values of {0{9t-)) m can be obtained. The general 
trend is around {O{9t-))o = 0, though. Then even for large M one must be careful with the 
computation of the integrated response function since a large initial error of x(0) due to the use 
of the average (. . .)m may significantly propagate as t increases. 

In order to avoid the problem of the sensitivity to the initial condition, instead of using 
directly the average (. . .)m in Eq. ([5|), one can define an estimator xa/(^ — i-u) satisfying the 
initial condition Xa/(0) = as required ideally by Eq. In this way the propagation of the 
initial error given by {0{9ti))M is suppressed at the beginning. An intuitive way to define xm 
can be outlined from the usual protocol to compute the integrated response function 

X{t-U) = jj (11) 

where the time t* is chosen such that 0{9f) = 0^^{9t-) and 0^^{9t) denotes the observable 
measured during the perturbed process. Notice that EqlTl]is justified by the fact that in the case 
of an infinite number of samples {O{9t))o = Vt, because of the time translational invariance of 
the NESS. In contrast when M is finite, it is useful to take into account, in EqlTll that (. . .)sa 
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and (. . .)o are performed independently the first on the perturbed trajectory 0^^{0{6t) and the 
second on the unperturbed ones O{0(t-ti+f)), specifically 



XM{t - U) 
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ti+f)) 



k=l 



where L is the number of unperturbed trajectories such Ok{ 
can be rewritten 



= 0' 



SAi 



j=l (. k=l ) 



(12) 



J. Therefore EqIT2] 
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where 5Ojk{0t) = O' 



turbed trajectory 0*''^( 



%) — Ok{0{t-ti+f)) is the the instantaneous difference between a per- 
\) and an unperturbed one Ok{0t). An example of this procedure is 
depicted in Fig. |31[a), where ti has been set equal to zero. We see that for a given perturbed 
trajectory O^"*(0t) (thick dashed red line), obtained after that 5A has been applied, one should 



look for an unperturbed NESS trajectory Ok{0{t)) such that Ok{0f) 
unperturbed trajectories shown by the solid lines. 
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Figure 3: (a) Examples of perturbed (thick dashed red line) and unperturbed (solid lines) 
trajectories used to estimate xi^) using Eq. (jl3p . Inset: expanded view at short time, (b) 
Estimate of -x{t) for M = 50, 100, 250, 500 and L = 1 (sohd hues) and M = 500, L = 200 
(dashed line). Inset: Comparison of xa/(0 for M = 500, L = 200 with the poor estimate done 
by the uncorrected average of Eq. ([T0| for M = 500. (c) Comparison between the experimental 
C(0) — C{t) and the best estimate of —[x{t) ± o'x(^)]- 

In this way 5Ojk{0ti) = by construction and the estimator defined by Eq. (|13p satisfies the 
condition xm{0) — 0. For M,L oo, x&iit — ti) converges to x defined by Eq. ([5]), because 
{{0{9t)))M — > 0. In Fig. [3lb) we show XM{t), with ti redefined as ti = 0, computed using 
Eq. (|I3p for different values of M and for fixed L = I (solid lines) and L = 200 (dashed line). As 
M increases for L ^ 1 the curves converge to a single profile which must correspond to that of 
x{t) ideally given by Eq. ([5]). The additional conditional average done for L = 200 smoothes the 
slightly fluctuating profile for M = 500 (thick solid blue line) resulting in the thick dashed solid 
line. For comparison we also show in the inset of Fig. [3jb) the raw estimate of x obtained using 
the average of Eq. ([T0| for the same M = 500 perturbed trajectories without correcting the 
effect of the initial sampling. In this case the propagation of the initial large error of {O{6o))o 
gives rise to a very poor estimate of the integrated response function for t > 0. 
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4.3 Experimental test 

Finally we proceed to test the theoretical fluctuation-dissipation relation ^ for the experimental 
unperturbed NESS trajectories of the Brownian particle and those perturbed around NESS. For 
this purpose we compare the best estimate XM{t) of x{t) done for M = 500 and Lj = 200 
with the right-hand side of Eq. The involved correlation function C{t) = {0{9t)0{9o))o on 
the right-hand side is computed using unperturbed NESS trajectories. In Fig. [3l^c) we compare 
XM{t) with C{t) — C(0). Besides, one can estimate the statistical error of the experimental 
X{t) at each f > by computing the standard deviation of {0^^{6t) — 0{9t))M over the Lj 
possible choices of the unperturbed Ok (t) . The standard deviation ±ay, (i) obtained in this way 
is also shown in Fig. [3]Jc) showing that after following the careful procedure to estimate x the 
relation x{t) ~ C!{t) — C(0) is verified by the finite experimental data. Note that without the 
finite correction {{0{6t))) m of the initial condition in Eq. one would largely underestimate 
the direct measurement of the integrated response function leading to an apparent violation of 
Eq. dH]). The results of the present paper are consistent with those of Refs. [T3J [TB] where two 
fiuctuation-dissipation formulae equivalent to Eq. ^ but involving different observables from 
the one studied here are checked experimentally. 

5 Conclusion 

We have studied the influence of finite sampling in the computation of the linear response 
function of a Brownian particle in a toroidal optical trap around a NESS. We have shown that 
there are two different effects that may lead to a very poor estimate of the experimental linear 
response function when the data analysis is not performed carefully. This is an important point 
that must be assessed in general when applying in experiment and numerical simulations the 
different generalized fluctuation-dissipation formulae recently derived for NESS. 
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